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I. INTRODUCTION AND OVERVIEW 
•% 

System identification technology has been used success- 
fully for many vehicles. Because of their large number of 
degrees of freedom and complex aerodynamic interactions, 
the rotorcraft have always presented a special challenge to 
system identification methods. A completely integrated 
methodology has been developed under this NASA contract to solve 
this difficult problem. This methodology has also been 
translated into a user oriented series of computer programs. 

This volume provides basic guidelines for efficient and 
effective use of one of these computer programs. 

Figure 1 shows a schematic flowchart of the overall 
data processing technique for rotorcraft. The first step in 
this procedure is state estimation and instrument calibration. 
This is implemented by the computer program DEKFIS (for 
Discrete Extended Kalman Filter and Smoother) which implements 
an extended Kalman filter/smoother using the Friedland -Duffy 
formulation. Instrument biases and scale factors are estimated 
at this stage together with any state which is not measured 
directly. The second step involves estimation of the mathe- 
matical model of various forces, moments and interchanges. 

This is implemented in OSR (Optimal Subset Regression) computer 
program which uses a regression technique. Accurate estimates 
of parameters are obtained in the final step. One of two 
computer programs is used for this . purpose . SCIDNT implements 
the maximum likelihood method for linear systems and NLSCIDNT 
extends the method to nonlinear rotorcraft models. 


The contract research effort which led to the results in this report was 
financially supported by the Structures Laboratory, USARTL (AVRADCOM), 
NASA Langley Research Center and NASA Ames Research Center. 
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Figure 1 Integrated Rotorcraft System 
Identification Procedure 


Accuracy of parameter estimates may be improved by using 
flight test inputs based on the input design program, INDES. 

This user's manual describes the SCIDNT computer program. 
The details of the theory and the particular implementation 
used are given in the final report.* 


Hall, W.E., Gupta, N.K., Hansen, R. and Bohn, J. , "State 
Estimation and Parameter Identification for Rotorcraft," 
Final Report on Contract NASI- 14549, May 1978. 
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I I . BACKGROUND 


The computer program Linear SCIDNT (from SCI maximum 
likelihood IDenTif ication for linear systems) evaluates rotor- 
craft stability and control coefficients from flight or wind 
tunnel test data. It implements the maximum likelihood method 
to maximize the likelihood function of the parameters based on 
measured input/output time histories. Linear SCIDNT may be 
applied to systems modeled by linear constant-coefficient 
differential equations. This restriction in scope allows the 
application of several analytical results which simplify the 
computation and improve its efficiency over the general non- 
linear case. The functions of Linear SCIDNT may be summarized 
as follows. For the linear system: 


^ x (t ) = F (e)x(t) ■+ G(0)u(t) + rce)w(t) 

y(t k ) = H(e)x(t k ) + D(9)u(t k ) + v.(t k ) ; k- 1 , 2 , . ... ,n 

where 

E [w(t)w T (t) ] - Q (9 ) <S (t-t) 

E[v(t k )v T (t jl )] = R(9) $ k)l 
E [v (t)w T (x) ] = 0 

9 = a vector of parameters to be identified 


where 

SCIDNT 


<5 (t-x) and 5 


kJl 


are Dirac delta functions. Linear 


estimates specified elements of the parameter vector 
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The program has the following special features: 

(1) The program is set up for rotorcraft models. 

(2) The applicable model may be specified depending on the 
measured response (this includes rotor, body, rotor/ 
body, coupled and uncoupled, longitudinal and lateral 
models) . The set of measurements may also be selected 
depending on available instruments . 

(3) The program estimates the standard deviations and con- 
fidence levels of all parameter estimates. 

(4) Up to 120 parameters can be identified in one run. 

(5) Multiple maneuvers can be processed in one run. 

(6) The user can specify that any of the parameters are 
known constants or are to be identified. 

(7) Extensive diagnostic printouts can be switched on to 
aid the user in setting up the input deck properly 
before making a complete run. 

(8) The user has the option of producing time history 
records of actual and estimated measurements in the 
form of tabular printout, printer plots, or magnetic 
tape for off-line high resolution plotting. 
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III. PROGRAM STRUCTURE 


This section gives a brief discussion of the major sub- 
routines and functional blocks of Linear SCIDNT. Figure 3.1 
shows the basic program structure. 

The main routine, DRIVER, is responsible for setting up 
the problem for the optimization routine SMAIN, calling SMAIN 
and then performing output tasks such as printing and/or storing 
(for later plotting) control and actual vs. estimated time 
histories. 

SMAIN, the optimization routine, uses the Levenberg- 
Marquardt search procedure to drive the identified parameters 
to values such that the likelihood function is maximized. As 
the parameters are stepped, information is printed out describ- 
ing the progress of each iteration. 

UPDATE is called by SMAIN to calculate the likelihood func- 
tion, its gradient, with respect to the parameters and the 
information matrix (an approximation to the Hessian of the like- 
lihood function with respect to the parameters). The cost and 
gradient are derived from propagation of the measurement and 
measurement sensitivity (with respect to each parameter) equa- 
tions, performed by UPDATE, using the transition matrix tech- 
nique. UPDATE also contains special code for the case when the 
model is assumed to have process noise excitation. 

In addition to the functional blocks described above, 
numerous other small routines are called to perform utility 
tasks, such as specially formatted printout, pointer manipu- 
lation, and matrix algebra. 


5 
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Figure 3.1 General Flowchart of SCIDNT 



IV. INPUTS TO THE PROGRAM 


Linear SCIDNT requires two classes of input. The first 
type, which will be referred to as "card- input" defines the 
model, its size (number of states, etc.), and its parameters, 
including indicators as to which ones are to be identified. 

Other quantities read include controls over the optimization 
routine, printout flags, and plot control flags. Card-input 
is read from unit 5. 

The second type of input referred to as "time-history 
input" consists of tabular values of the measurement and control 
time histories. SCIDNT calls a subroutine INREAD once before 
beginning the identification algorithm. Subroutine INREAD 
reads the values of the measurements y and controls u for 
the entire time period of the experiment. This is because in 
general, INREAD must read data in many different formats for 
various types of simulation or flight test data. 

4.1 CARD INPUT 

4 . 1 .1 General 

Card-input to Linear SCIDNT consists of three different 
groups of cards. The first group, consisting of only one card, 
specifies a 76- Character title to be printed at the top of 
every page of output. This serves to identify the run and/or 
the data which is specified as input. 

The second group, a NAMELIST data deck called INPUT, speci- 
fies the problem size (number of states, number of measurements, 
etc.), optimization algorithm controls, and program output 
controls. 
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The third group of cards specifies the numerical values 
of the parameters in the model and flags the ones to be iden- 
tified. Parameters not flagged are assumed to be fixed and 
known. Also indicated on these cards are parameter labels, 
parameter initial values, and upper and lower bounds for iden- 
tified parameters. The bounds are useful for parameters such 
as sines or cosines, whose magnitude must not exceed 1. 

A detailed description of card input is seen in Table 4.1. 
The NAMELIST cards are punched according to usual NAMELIST 
conventions, such as free format between columns 2 and 72, order 
independence, etc. See the applicable FORTRAN manual for more 
details. A sample input deck is shown in Figure 4.1. 

4.1.2 Detailed Explanation of Selected Card-Inputs 

K1MAX 

K1MAX specifies the maximum number of successful itera- 
tions that the optimization algorithm will perform before 
returning to the main program. Each successful iteration 
represents a step toward the maximum of the likelihood function 
(equivalently, the minimum of the negative log of the likelihood 
function, which is the actual cost function used) . If conver- 
gence is achieved in less than K1MAX iterations, then the final 
values of the parameters will reflect the last step calculated. 
If convergence is not achieved in K1MAX iterations, the para- 
meters will have the values before the last step was calculated. 
This feature allows the initial and final parameter values to 
be equal when K1MAX is input as 0. In this fashion, time 
history plots of the system with the initial values of the para- 
meters can be generated, and the first step calculated on a low- 
cost checkout run. An even lower-cost checkout run can be made 
by setting K1MAX = - 1 , which is the same as setting KlMAX * 0 , 
but without a calculation of the first step. 



Table 4.1 
Input Card Formats 


CARD GROUP 

COLUMNS 

FORMAT 

VARIABLE NAME AND DESCRIPTION 

Title 

1-76 

19A4 

Run and/or Program identification. 

INPUT 
Name! ist 

2-72 

Free 

First card must begin with $INPUT followed 
by SEND. Intervening cards set values for 
the following variables (defaults in 
parentheses) r 



Integer 

NS - number of states (22) 



Integer 

NP 3 number of measurements (25) 



Integer 

NQ 3 number of controls (8) 



Integer 

NN 3 number of time points to be read by 
INREAD from data file (21). Can be 
negated to trigger data printout by 
INREAD. 



Integer 

K1MAX 3 maximum number of iterations (6). 
See notes. 



Integer 

K2MAX * maximum number of cost function 
increases per iteration (4). See 
4.1.2. (step cuts) 

■ ■ 


Integer 

MCYCLE 3 number of parameters identified 
per iteration (min(# parameters 
i d ' d » 12 ) ) . Must be less than or 
equal to total number of para- 
meters identified. See 4.1.2. 



Integer 

NG = number of process noise sources (6). 
If equal to zero, output error 
option is assumed. See A. 1.2. 

' 


Integer 

IRCMP = 0, R-matrix not updated. 

IRCMP > 0, R-matrix updated every IRCMP . 
iterations. See 4.1.2. 




Table 4.1 (Continued) 


CARD GROUP COLUMNS FORMAT VARIABLE NAME AND DESCRIPTION 

INPUT 2-72 Free 

Namelist 

(Cont'd) Integer IPRT = 0, No diagnostic printout 

IPRT = 1, State and measurement debugging 
printout 

IPRT = 2, State and measurement sensitivity 
debugging printout 

IPRT = 4, Information matrix debugging 
printout. 

Default = 0. 

Note: These values may be summed to obtain 

combined printouts. Also, the sum 
may be negated to obtain time his- 
tories associated with each level 
of printout. See 4.1.2. 


Real DELTA = Sample time increment of input 

measurement and control time his- 
tories (.05 seconds) 

Logical EIGF = Calculate and print out F-matrix 
eigenvalues before and after 
optimization procedure (.FALSE.) 

Logical PRTTB = Print out tabular time history of 
controls and actual vs. estimated 
measurements. (.FALSE.) 

Logical PRTPL = Produce printer plots of controls 

and actual vs. estimated measure- 
ments (.FALSE.) 

Logical TAPEPL = Produce mass storage file of con- 
trols and actual vs. estimated 
measurements. These time histor- 
ies are written on unit 2 and can 
be used as input for off-line 
plotting. (.FALSE.) 

Integer NREC = Number of maneuvers in data record 

0 ) 

Integer IREC = The Ith maneuver starts at the 
Arra y IREC(I)th data point (1) 






Table 4 . 



COLUMNS 

INPUT 2-72 Free 

Namelist 

(Cont'd) Integer 

Integer 

Integer 


CARO GROUP 


FORMAT 


Logical 

Real 

Real 

Real 

Real 

Real 

Real 

USER LABELS 1 A1 

3-4 12 

11-50 4A10 


(Continued) 


VARIABLE NAME AND DESCRIPTION 


MAXP * maximum number of parameters 
(computed from other inputs) 

NTERMS = # of terms in transition matrix 
series evaluation (20) 

IPRTSM * 0 prints inputs for FMARO, PMARQ, 
NTERMS, odd turns on level 1 
print (0) 

USERLAB = .TRUE, to read user supplied 
labels for plots (.FALSE.) 


FMARQ = initial Marquard search 
parameter (30. ) 

PMARO = factor for increase or decrease 
of Marquard parameter ( 1 . ) 

DEFBDU 35 upper bound default on P (10 6 ) 

DEFBDL = lower bound default on P (-10®) 


STEPMN > 



def i nes; convergence 


RELERR > | relative change in likelihood 

functionj defines convergence (10" b ) 


CHAR * Y for measurement label 
* U for control label 

INDEX = which measurement or control 

ALABEL (I), I = 1,4 label 
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Table 4. 


CARD GROUP 

COLUMNS 

FORMAT 

Parameter 

1 

A1 


2 

A1 


6-25 

415 


26-35 

A10 

. 

36-50 

E15.7 


51-65 

E15.7 


66-80 

E15.7 


(Continued) 


VARIABLE NAME AND DESCRIPTION 


J2 » blank: this parameter is to be held 

fixed throughout the optimi- 
zation 

J2 * non-blank: this parameter is to be 

identified 

ECHK = non-blank: terminator card for the 

parameter group 

(IPX(I) ,1=1,4), pointers to positions in 
which this parameter appears 
in linear system (see 4. 1.2) 

PLAB: Parameter name. 

PVAL: If this parameter is identified, P 

is its starting value. If it is 
not identified, P is its known 
value. (Default = 0.) Fixed para- 
meters of value zero need not be 
specified in the parameter card 
group. 

PLJ1: Lower bound of this parameter, if 

identified. Default = -1.0E+6 

PUJ1: Upper bound of this parameter, if 

identified. Default = 1.0E+6 


See 4.1.2 for more detailed description 
of parameter card group 
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COLUMN 

1 


ROTORCRAFT SCIONT TEST CASE 

* INPUT 

A *': 

NP«l*. 

N <) a *. 

' 10*0 • 

«* APUa 30.. P » APO » l ., 
NIP** * 10. 

K2H**»?0. 

NN »501 » 

DELTA*. 0] . 
TAPEPL».TPUE.. 
EIOFa.TRUF.. 
P»TPL*.TPUE.. 

SENO 

253 
25* 

255 

256 

257 
25 * 


1 

2 

3 

4 

5 

6 
10 
12 
1* 

15 

16 
17 
16 

19 

20 


267 

26 * 

269 

270 

271 

272 


Title Card 


SINPUT Namelist 


-.00*6 
.0156 
•. 12 * 
-.00 11 
-.0006 
-.0016 
.029* 
. 0*91 

•,0a** 

. 000*2 

-.0667 

-.0071 

-.015* 

.0023 

.0061 


PARAMETER SPECIFICATION 
CAROS 


10 Flag 


2* 


BO V 

.0005 

26 


9 A V 

-.0707 

26 


*HV 

-0.0506 

29 

2*1 

XV.’ 

.09*2 

30 

2*2 

Y* 

.00216 

31 

2*3 

z* 

9,7*. 

32 

2** 

t> / 

-.0099 

33 

2*5 


.0016 

3* 

2*6 

NW 

-.0002 

3* 


«0» 

.*32 

*0 


BA* 

.0379 

*2 


HeU 

-.09*5 

*3 

295 

XP 

-.126 

** 


YP-VSALO 

6.63 

*5 

297 

ZP. 

-6.12 

*6 

29* 

LP 

-.172 

*7 

299 

MB' 

-.0055 

k* 

300 

NP 

.111 

*6 

135 1*6 195 

ON£ 

1. 

5? 


eiOP 

.*96 

5* 


PAP 

—6 .*6 

5.6 



-16.2 

57 


xo-vstho 

-9.9 7 

6* 

310 

YO 

-.P69 

a9 


za*vCTM0 

193. 

*0 

312 

LG 

.017 

bl 

313 

*<i 

-5.0 

6? 

31* 

Mi 

.069* 

63 


SPHOTTHO 

0. 

6* 


CPHO 

1. 

66 


BOO 

-3.97 

4* 


KAU 

— 1 * . * 

•arameter Location 
Pol nters 

Parameter 

Name 

Startup Values 


Figure 


4.1 Sample Input Data 
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70 



12.1 

71 

323 

X8 


7? 


YP-VCAIO 

-lo7.0 

73 

325 

70 

*.* 

7* 

326 

L* 

.**3 

7* 

327 

M8 

.00*5 

7* 

328 

NH 

-.726 

77 


CP80TTM0 

.0**7 

78 


-SPM0 

0. 

aO 


808 

-3.96 

HP 


9*8 

.156 

8* 


889 

-.72 

88 


OCPMOC7MO 

J2.P 

87 


-8SPH0CTH0 

0. 

99 


-6CTH0 

-32.2 

100 


-9SPH0ST80 

0. 

101 


-«CPHOSTMO 

-1 .»* 

113 

365 

XAO 

25.1 

114 

366 

7*0 

•10. 1 

US 

367 

ZAO 

-21a. 

118 

368 

L*0 

-.52 

117 

369 

9*0 

-.06* 

U« 

370 

ftao 

-l.» 

122 


80*0 

-166. 

12* 


8**0 

89.6 

126 


88*0 

-8.31 

127 

379 

XAOOOT 

• 0* 1 6 

128 

380 

Y AOOOT 

-.376 

129 

381 

ZAOOOT 

.799 

130 

382 

laooot 

-.01 

131 

383 

M A 000 7 

-.0-21 

132 

38* 

NAOOOT 

-.0152 

136 


80*0007 

-13.7 

138 


gAAOOOT 

3.32 

1*0 


B8A00O7 

1.56 

1*1 

393 

X* l 

-2*. 

1*2 

39 * 

Y A 1 

-.3? 

1*3 

389 

Z*l 

-11. 1 

1** 

396 

.LAI 

-*.56 

1*9 

397 

MAI 

3.8* 

1*6 

398 

NA1 

-.83 

150 


50*1 

19.0 

152 


8**1 

-*7. 

15* 


88A1 

15*. 

155 

*07 

XA100T 

-2.51 

156 

*08 

YAIDOT 

-.21 

157 

*09 

Z 4 1007 

3.78 

158 

*10 

L»1D0T 

-.*72 

155 

*11 

MA 1007 

.3 

160 

*12 

6A100T 

-.0532 

16* 


80A1DOT 

.999 

166 


PAAIOOT 

-19.* 

166 


an* 1007 

22.8 

16* 

*21 

XSl 

-1.91 

170 

*22 

yPi 

25.8 

171 

*23 

Z81 

11.7 

17? 

*2* 

L«1 

17.8 

173 

*25 

MH 1 

1.0 

17* 

*26 

NPl 

1.8* 

178 


a 081 

-I6.e 

180 


8Aai 

•161. 

182 


889 1 

-48 . 7 

183 

*35 

XMIOOT 

-.45 

18* 

*36 

YhlOOT 

.923 

186 

*37 

ZH10OT 

-*.»* 

186 

*38 

LilQOT 

1 .07 


parameter 

SPECIFICATION 
CAROS (CONT.) 


Figure 4.1 (Continued) 



187 

*39 


««tOOT 

.10* PARAMETER 

188 

**0 


MritDOT 

.08 SPECIFICATION 

19? 



■ aoaiooT 

-•59 CAROS (CONT) 

19* 



9*9 10OT 

-23.7 

198 



968100T 

-17.6 

296 



yP 

-.91 

309 



XQ 

-2. *2 

197 

4*9 


X TH IS 

-*.92 

198 

*50 


tTHIS 

12.0 

199 

*51 


2ThIS 

12*. 

200 

*52 


ltkis 

.913 

?01 

*53 


-This 

-.352 

?0? 

*5* 


NThlS 

-.0905 

208 



80TH1S 

— 6* . 2 

208 



84TH1S 

-117. 

210 



M8TH1S 

163. 

211 

*83 


XTHIC 

-*.92 

?1 2 

*6* 


YTWIC 

11.9 

213 

*65 


2TH1C 

17.1 

?V* 

*66 


ltwic 

3.33 

215 

*67 


mTMIC 

-.0796 

218 

*68 


MH1C 

-.5*3 

220 



90TKIC 

9.82 

222 



RATH1C 

163. 

22* 



88T81C 

110. 

225 

*77 


XOLTR 

-.*8 

229 

*78 


yoltr 

6.«7 

227 

*79 


ZOLTR 

.53* 

228 

*80 


LOLT9 

3.65 

229 

*81 


uOLTR 

.0078 

230 

*82 


nOLTR 

. -5.17 

23* 



RODLTB 

-1.9* 

238 



RAOLTR 

-.63 

238 



98ULT9 

.29 

239 

•91 


XTWO 

26. S 

2*0 

*92 


YTHO 

-.706 

2*1 

*93 


ZT80 

-269. 

2*2 

*9* 


ltho 

-3. *5 

2*3 

*95 


mTMO 

,66 

2** 

*96 


MT>0 . 

*.7 . 

2*8 



ROThO 

231. 

250 



B«TH0 

30.7 

252 



BHTHQ 

-8*. 

311 



Z9 

2*.* 

32* 



Y8 

1 .61 

301 

316 

331 

3*6 ON* 

l. 

381 

376 

*05 

*3* ON6 

1. 

*09 



WAX 

1.256-02 

920 



SAY 

1.256-02 

*35 



B»Z 

1.256-02 

550 



UPDOT 

1.636-05 

595 



PUOOT 

1.536-05 

580 



UROOT 

1.536-05 

595 



SP 

3.796-06 

810 



R9 

3.796-06 

825 



UR 

3.796-06 

6*0 



RP>4! 

1. *56-07 

655 



RT86TA 

1.456-07 

670 



RAO 

5.786-07 

685 



9A 1 

5.786-07 

700 



09 1 

5.766-07 


Figure 4.1 (Continued) 



K2MAX 


Sometimes the calculated step of a given iteration is too 
large, causing the negative log likelihood function to increase. 
When this happens, the optimization algorithm cuts back the 
step until the cost decreases. K2MAX limits the number of 
times that this can occur in an iteration. The default value 
of 4 should be sufficient for most cases. 

MCYCLE 

Although every parameter can be moved during each itera- 
tion, significant savings of computer time can be realized if 
only a portion of the total number of identified parameters are 
stepped each iteration, especially if there are many (> 20) 
parameters to be identified. MCYCLE determines how many para- 
meters are stepped at each iteration. If MCYCLE =0, or is 
absent, then all the parameters are stepped at once each itera- 
tion, up to a maximum of 60. If MCYCLE > 0, then MCYCLE para- 
meters are stepped. The identified parameters are grouped 
MCYCLE at a time in the order that they appear in the parameter 
cards. Different MCYCLE groupings can be achieved by re- 
arranging the parameter cards. If the total number of identi- 
fied parameters is not an even multiple of MCYCLE, the remainder 
of parameters in the last MCYCLE group are taken from the first 
group, and subsequent MCYCLE groups are changed accordingly. 
MCYCLE must be less than or equal to the number of identified 
parameters, and must never be greater than 60. 

Note that this advanced computation technique was intro- 
duced into the computer program to save execution time while 
solving large order problems, the final estimates of para- 
meters should be the same (within limits of numerical error 
and one standard deviation of parameter estimation error) , 
irrespective of the value of MCYCLE. This technique has been 
used in the optimization of large static econometric models 
and is based on a similar technique proposed by Golub and 
Pereyra [1] . 
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NG 

Setting NG > 0 activates the process noise option of Linear 
SCIDNT, in which the linear model is assumed to be excited by 
random disturbances. If this option is used, parameters for 
r and Q matrices should appear in the parameter cards. 

Since the process noise case assumes that the estimated state 
is the output of a Kalman filter, which must be calculated, 
significant computational overhead is required for this option. 
However, improved accuracy of parameter estimates can be 
realized for data contaminated with process noise. If NG > 0, 
the input values of the R and Q matrices determine the 
dynamics of the state estimator. 

If NG is set to 0, the output error option of the program 
is used. In this case, all errors between estimated and 
observed measurements are assumed to be due to measurement 
noise. No T or Q matrices are used, and the input R 
matrix is not used in the calculations, as measurement noise 
covariance is estimated from the data. Output, error mode runs 
faster than process noise, and should be used when disturbances 
are known to have negligible effect on the data. 

IRCMP 

In the process noise option, the input value of the R- 
matrix is used to compute the Kalman filter gains. However, if 
the program has run for a few iterations, a better estimate of 
R can sometimes be obtained from the data and from new esti- 
mates of the process noise covariance. If IRCMP > 0, then the 
R-matrix used in the Kalman filter calculation is updated by 
the program every IRCMP iterations. If IRCMP = 0, R is held 
fixed at its input value. Every time R is updated, its new 
value is printed out. 

In the output error option, IRCMP has no effect on the 
calculations. However, better estimates of R for future runs 
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can be obtained by setting IRCMP « 1, which will cause an esti- 
mate of R to be printed at each iteration. 

IPRT 

IPRT can be used to trigger several levels of diagnostic 
printout. Positive values of IPRT display matrices associated 
with each level of printout. (Level 4 is the same as Level 2, 
since there are no new matrices associated with the information 
matrix calculation.) Level numbers can be summed to obtain 
combined printouts. For example, IPRT * 3 causes a combined 
Level 1 and Level 2 printout, and IPRT * 7 includes all levels 
of printout. If the value of IPRT is negated, then time 
histories associated with each level are printed, as well as 
matrices. (IPRT * -4, -5, -6, or -7 will include a time 
history printout of the sensitivities used to compute each 
element of the information matrix) . 

Caution : Values of IPRT of magnitude greater than 1 will 

cause voluminous output. If such printouts are desired, it is 
best to set K1MAX » 0, and identify just a few parameters 
(1 or 2, if possible). If IPRT < 0, try to limit the number of 
data points to be 20 or fewer. 

IPX Vector 

IPX is read from a parameter card and determines the posi- 
tion in which the parameter appears in the system matrices. 

A parameter may appear in as many as four matrix elements, 
and the pointers to these elements are determined as follows. 

If the columns of the matrices (in the order F, G, H, D, R, r, 

Q) are stacked into one long vector, then the pointer to a 
given matrix element is equal to that element's position in 
the long vector. The four elements of IPX correspond to the 
four matrix elements in which a parameter can appear. Of the 
four fields on the parameter card for IPX elements, the leftmost 
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ones are used for non-zero values of IPX, i.e., unused IPX 
fields must appear on the right. The ranges of values for the 
pointers are as follows: 

Elements 


In 

Go From 

TO 

F 

1 (IF) 

IF+NS*NS = (IG- 1) 

G 

IG 

IG+NS*NQ = (IH- 1) 

H 

IH 

IH+NP*NS = (ID-1) 

D 

ID 

ID+NP*NQ = (ID-1) 

R 

IR 

IR+NP*NQ -(IGAM-1) 

r* 

IGAM 

IGAM+NS*NG = (IQ-1) 

Q* 

IQ 

IQ+NG*NG 


4.2 TIME HISTORY INPUT 

Time history input is read by subroutine INREAD. Because 
measurement and control time histories can be retrieved from 
tape, cards, or disk in numerous formats, INREAD is left as a 
user-written routine. The INREAD used in the sample run is 
shown in Figure 4.2. 

Subroutine INREAD 

SUBROUTINE INREAD (NP ,NQ ,NN, DELTA) 

COMMON/TLCM/T (1) 

COMMON/ULCM/U (1) 

COMMON/YLCM/Y (1 ) 

LEVEL 2, T, U, Y 

Inputs: NP = number of measurements 

NQ = number of controls 

NN * number of data points 

DELTA = time history sample interval (seconds) 


*Used only with process noise option 
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Outputs: 


Notes : 


T * array for storing the time of each data 
point 

U * array for storing the control time history 

Y * array for storing the measurement time 
history 

NN = number of data points actually read from 
file 

(1) The program uses the array T only for the plot 
outputs. The identification algorithm does 
not require T. The variable DELTA is passed 

to INREAD mainly for calculation of this array, 
and otherwise may be ignored (although a 
place must be held for it in the parameter 
list). This calculation can be done if the 
data file does not explicitly have on it the 
time for each data point. 

(2) As noted under "outputs," NN may be reset by 
INREAD to reflect the actual number of data 
points read from the file. This must be done 
if, for example, an end-of-file is reached 
before the anticipated number of data points 
are read. If NN is reset, a message to that 
effect should be printed. 

(3) INREAD should not reset the values of NP or 
NQ. DELTA may be reset, but if that occurs, 

a message to that effect should be printed, as 
an incorrect sample time can have disastrous 
effects on later computation. 

(4) INREAD may perform any desired pre-processing 
of Y and U. 

(5) It is recommended that INREAD print out the 
Y and U as arrays immediately after reading 
them to check the validity of the data, as shown 
in the example. DRIVER, which calls INREAD, 
has been coded so that a negative value of NN 
can be read from NAMELIST $INPUT, and then 
passed to INREAD to trigger this printout, 
should the user desire to incorporate it. 

If NN comes into INREAD as a negative value, 
it must be reset positive via the IABS func- 
tion before returning to DRIVER. 
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C6) If the user should require additional card 
input in INREAD to control the processing of 
the input data, the cards can be placed 
immediately after the parameter description 
group of cards. Card input is on unit 5. 

(7) The arrays T, U, and Y are in LCM. Note that 
the INREAD declaratives specify only the 
starting addresses of these arrays and not 
their dimensions. Therefore, when these 
arrays are read they must be considered singly 
dimensioned, and the variables NP, NQ, and NN 
must be used to calculate the proper sub- 
script for these arrays. T is a Vector, U is 
assumed to have NQ rows and NN columns, and 
Y is assumed to have NP rows and NN columns. 




V. OUTPUTS OF THE PROGRAM 


Four forms of output can be obtained from Linear SCIDNT. 
They are: 

(1) Tabular printout showing the progress of the parameter 
estimation process. 

(2) Tabular printout of time histories of controls, 
measurements, and measurement estimates using the 
final values of identified parameters , (PRTTB = .TRUE.) 

(3) Printer plots of control, measurement, and measure- 
ment estimate time histories, (PRTPL = .TRUE.) 

(4) Control, measurement, and estimated measurement time 
history data stored on mass storage to generate off- 
line plots, (TAPEPL » .TRUE.) 

The first type of output is always produced. The last three 
are optional and may be independently selected by the user. 


5.1 TABULAR PRINTOUT OF ITERATION PROGRESS 


Figure 5.1 is an example printout of a SCIDNT run using 
the card input shown in Figure 4.1. Explanations of the results 
are noted on the printout. 
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5. 2 PRINTER PLOTS 

Figures 5.2 and 5.3 show typical printer plots of measure- 
ment and control time histories. Actual and estimated time 
histories for each measurement are plotted on the same axis 
sd that closeness of fit can be visually evaluated. 

5. 3 TIME HISTORY ON MASS STORAGE 


This option is provided to put measurement and control 
tile history data on a mass storage file, primarily for the 
purpose of producing off-line plots analogous to the printer 
plots. Data is written to the device (logical unit 2) in 
unformatted records as follows: 


Record 


Data 


1 


2 through NN+1 


NN+2 


NS, NQ, NP, NN 

NS - number of states 

NQ * number of controls 

NP 3 number of measurements 

NN 3 numb e r of data points per variable 

T(J) , (U (K, J) ,K 3 1 ,NQ) , (Y (K > J) ,K*1,NP) , 
(YHAT (K, J) , K*1 ,NP) 

T(J) = time for Jth data point 

U (K, J) 3 Kth control at Jth data 
point 

Y (K, J) 3 Kth actual measurement at 
Jth data point 

YHAT (K, J) * Kth estimated measurement at 
Jth data point 

End-of-file mark 
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Figure 5.2 Measurement Time History 
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Figure 5.3 Control Time History 



APPENDIX A 


A. 1 MAXIMUM PROBLEM DIMENSIONS 

The Ames 7600 implementation of Linear SCIDNT is dimen- 
sioned to the following maximum problem sizes: 

NS (number of states) _< 22 

NQ (number of controls) _< 8 

NP (number of measurements) _< 25 

NG (number of process noise sources) <6 

Total number of identified parameters £ 120 

MCYCLE (number of parameters identified per iteration) _< 60 

NN (number of data points) <_ 1001 

A. 2 JOB CONTROL LANGUAGE REQUIRED TO USE LINEAR SCIDNT ON 
AMES 7600 

The following job deck was used to make sample runs on 
simulated data. It can be used as an example for other SCIDNT 
job setups. 

RCSCI ,Tm,Pn, YD1. Job card, m and n depend on the particu- 
lar run. 

ACCOUNT, XXXXXX,TXXXX. Account card. 

MOUNT, VSN=D0185A, FSDMOHA. Mount the private disk with the 

SCIDNT program and input data. 

SETNAME, FSDMOHA. Establish the private disk as the default 

setname for subsequent ATTACHes . 

ATTACH, RCCH5 3, ID=VERMONT. Attach the time history data 

file (card images in UPDATE format). 
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UPDATE,P»RCCHS3,C=TAPE9,D,F. Un-compress the data file. 
ATTACH, SCIDNT , LIDBN , ID* VERMONT. Attach the SCIDNT program. 
SCIDNT. Load and execute SCIDNT. 

7/8/9 

This null logical record is needed for the UPDATE 
processor. 

7/8/9 

SCIDNT data cards (see Figure 4.1) 

6/ 7/8/9 
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